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Abstract 

A complete derivation, from first principles, of the reaction-rate for- 
mula for a generic process taking place in a heat bath of finite volume is 
given. It is shown that the formula involves no finite-volume correction. 
Through perturbative diagrammatic analysis of the resultant formula, the 
detailed-balance formula is derived. The zero-temperature limit of the for- 
mula is discussed. Thermal cutting rules, which are introduced in previous 
work, are compared with those introduced by other authors. 
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1 Introduction 



Ultrarelativistic heavy-ion-collision experiments at CERN and RHIC lead us to enter- 
tain a hope of reviving quark-gluon plasma (QGP) in the present day. As promising 
observables of the QGP formation, rates of various reactions taking place in a QGP 
(heat bath) have been computed by many authors. Almost all of them, however, 
concentrated on the analyses of particle production from a QGP or the decay rate of 
a particle in a QGP, whose computational method has long been known . 

Since then, through analyses from first principles, a calculational scheme of the 
rate of a generic thermal reaction has been proposed |2], |3|, [|, |, || 0. The resul- 
tant reaction-rate formula is written in terms of the Keldish variant of the real-time 
formalism (RTF) of thermal field theory ||. However, complete analysis of classes 
of diagrams, which leads to diagrams in RTF including thermal propagators with n 
(> 2) thermal self-energy insertion, is still lacking. Ref. || is the only work that 
discusses such classes of diagrams in scalar field theory. In the course of deduction 
H of such diagrams, there comes about an involved series, for which an identity is 
assumed. As for [|J, where fermion fields are dealt with, the set of diagrams under 
consideration is not analyzed. This is also the case}] for 0. Incidentally, the ther- 
mal self-energy part in itself and the one thermal self-energy-inserted propagator are 
deduced in @, |, |,|]. 

The principal purpose of this paper is to present a complete derivation of the 
thermal reaction-rate formula (Sees. II - V). 

There has been confusion regarding the issue of finite-volume corrections to the 
standard thermal perturbation theory. (Why and how has the confusion arisen is 
described historically in ]TtJ with relevant references.) By employing a cubic system 
with periodic boundary condition, it has been shown in that thermal expectation 
values of normal-ordered products of field operators can be chosen to be zero and there 
is no finite-volume correction on thermal amplitudes. It should be stressed that this 
statement is the statement within the RTF. The statement does not tell us whether 
or not the thermal reaction-rate formula deduced from first principles is free from 
finite-volume corrections. We shall derive in Sees. II - V the thermal reaction-rate 



tin fact, in [J7|, an n (> 2) thermal self-energy-inserted propagator is not deduced from the 
starting formula but is assumed at the start to have the correct form in RTF (cf. Eq. (17) in 0). 
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formula for the finite-volume system and explicitly see that there is no finite-volume 
correction. 

It should be emphasized that the absence of finite-volume corrections here as well 
as in (JTOj] is of rather academic since a cubic system with periodic boundary condition 



is taken. For physical finite-volume system, there are [Q two sources of entering the 
finite-volume effects on the thermal perturbation theory constructed on the basis of 
(grand) canonical ensemble. The one comes from the physically sensible boundary 
condition on the single-particle wave function. The other comes from taking the 
physically sensible ensemble. For the case of nonequilibrium case, such as expanding 
QGP, the situation is of course much more involved. 

In Sec. VI, through diagrammatic analysis for the reaction-rate formula, we derive 
the detailed-balance formula. In Sec. VII, we analyze the zero-temperature limit of 
the reaction-rate formula and reproduce a variant of the Cutkosky rules . 

At zero temperature, the cutting (Cutkosky) rules [|TTJ are the powerful device to 
investigate the imaginary or absorptive part of a scattering amplitude and a reaction 
rate like a scattering cross section. Then, it is natural to infer that a finite-temperature 
extensions of the cutting rules (thermal cutting rules) also plays an important role in 
thermal field theory. 

Previously, several authors |2|, [| |, ||, ^ 0, ||, [13], [H], [15], [HJ have discussed 
thermal cutting rules.[] However, because of the fact that the generalization of the 
notion of "cutting" in vacuum theory to the case of thermal field theory is not unique, 
the terms "cutting" and "(un)cuttable" are endowed with different meanings in |~2]. 0. 
0, f|, |, |6|, |7|, |[| [L3|, [L4|, [15], which causes recent controversy. With this circumstances 
in mind, we pigeonhole different definitions of thermal cutting rules (Sec. VIII). 



2 Preliminary 

We consider a heat-bath system of temperature T, composed of the fields <jr a \ with 

a labeling collectively a field type and its internal degree of freedom. We assume 

T » m and ignore m (hot plasma). The system is inside a cube with volume 

* Relationship between a thermal self-energy part (in imaginary-time formalism) and a rate of 
decay (production) of a particle in (from) a heat bath was clarified in [fij, from which the cutting 
rules as applied to the self-energy part can be read off. 
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V = L 3 . Employing the periodic boundary conditions, we label the single-particle 
basis by its momentum pk = 2irk/L, kj = 0, ±1, ±2, ■ • • , ±oo (j = x, y, z). 
Physically interesting thermal reactions are of the following generic type, 

{A} + heat bath {B} + anything . (2.1) 

Here {A} and {B} designate group of particles, which are not thermalized, such as 
virtual photons and leptons. (Generalization to more general process, where among 
{A} and/or {B} are is straightforward and will be dealt with in Sec. V.) The 

reaction rate 1Z of the thermal process ( |2.1|) is expressed @, ||, (| as an statistical 
average of the transition probability W = S*S (with S the ^-matrix element) of the 
zero-temperature (T = 0) process, 

{A} + {nl a) }^ {i?} + {4 Q)/ }, (2-2) 

where {4°^} denotes the group of <f)( a ^s, which consists of the number nj^ of 
(0( Q ) in a mode k ): 

K = AT/V, (2.3) 
Kf - V ^(Process (2.2)) 

V = E {n (") } PE {n (")' } ^o({4 Q) }- W a) '}), (2-5) 
p = iV^expf-^^^n^pk) . (2.6) 

V a k / 

Here (3 = 1/T, p k = |Pk|, and 2tt5(0) — tf — tj (~ oo) is the time interval during 
which the interaction acts. Wo = SqSq, the "thermal vacuum bubble," is the T = 
transition probability of the process indicated in Eq. (|2.5|) , i.e., the reaction among 
the heat-bath particles ^"^'s alone. Note that the perturbation series for T> starts 
from 1, 

£> = ! + .... (2.7) 



In Eq. (|2.6|) , iV is the normalization factor. In Eqs. (|2.4|) and (|2.5| ), stands for 
the summation with symmetry factors being respected, and, for a bosonic (fermionic) 
4>( a \ runs over 0, 1, 2, ■ ■ • , oo (0 and 1). It is to be noted that {A} and {B} in S, 
which we write {A, B} s , are not necessarily involved in one connected part of 5*. This 
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is also the case for {A, B}s*. We assume that, in W = S*S, {A, B}s and {A, B}$* 
are involved in one connected part, which we simply refer to as connected W. Then, a 
connected W consists, in general, of two mutually disconnected parts, the one includes 
{A, B}s and {A, B}$* and the other is a group of spectator particles. Generalization 



to other cases is straightforward. Examples of double-cut diagrams [17] for S*S are 
depicted in Fig. 1. It should be remarked on the form of p in Eq. ( |2.6p . Let us recall 
the following two facts. On the one hand, the statistical ensemble is defined by the 
density matrix at the very initial time t{ (~ — oo). On the other hand, in constructing 
perturbative RTF, an adiabatic switching off of the interaction is required ][H|, |, |J. 
Then, the Hamiltonian H in p = N~ l e~P H should be the free Hamiltonian H , which 
leads to Eq. (£§). 

As will be seen below, diagrammatic analysis shows that Af, Eq. ( |2.4|) , takes the 
form, 

M = M con V, (2.8) 
where Af con corresponds to a connected diagram and V is as in Eq. (|2.$). Then 

7^- A/con* 

The T = 5-matrix element is obtained through an application of the reduction 
formula. As an illustration, we take a heat-bath system of thermal neutral scalars 
0's, and we take {A} to be {$(pj);z = 1, ...,m} and {B} to be {$(q.,);j = l,...,n} 
with <3> a nonthermalized heavy neutral scalar. Assuming a $-0 coupling to be of the 
form $0', we have || § 



m n 

s = n (^,0 n (^,0 n 



3=1 



k=l 



"k 



i k =0 i^=0 



^sr' ri (*^,„o ft (<^k,n) (o i t 



n'=l 



n=l 



n 

n'=l 



where 



In Eq. 




n*»n*iii** 

n=l j=l fc=l 



1/2 



0) 



4! ^ 



(2.9) 



(2.10) 



denotes the Kronecker's 5-symbol, 
1 



k,n ' ' ' 



2 Pk yz 4> 



d xe 



□ 
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2E j VZ< s> 

x(a + M 2 ) ■■■$ j (x), (2.11) 



where Ej = Jp 2 + M 2 with M the mass of $. Z's in Eq. (|2.11 ) are the wave- 



function renormalization constants, So in Wq = SqSq is given by a similar expression 
to Eq. Q2.9D , where factors related to the $ fields are deleted. It is to be noted that, in 
Eq. (|2.9|), among (n k ) of 0k's in the initial (final) state, «k (i k ) of k 's are absorbed 
in (emitted from) the (i' k ) vertices in S. Remaining n k — iy_ (= n k — i k ) of 0k's 
are merely spectators, which reflects only on the statistical factor in A in Eq. ( [3.1 4| ) 
below. 

The expression for S*, the complex conjugate of S, is obtained by taking the 
complex conjugate of Eq. (|2.9|) , where we make the substitution, 

ik -> A 4 Jk • 
This applies also to the expression for Sq. 

3 Derivation of the reaction-rate formula 

In this section, we take self-interacting neutral scalar theory. Generalization to the 
complex-scalar theory is straightforward (cf. Sec. VIII). A comment on gauge theories 
is made at the end of this section. Fermion fields are dealt with in Sec. IV. 

3.1 Analysis of non mode-overlapping diagrams, «k + *k + jk + 

Jk<2 

In this subsection, for completeness, we briefly recapitulate the heart of the analysis 
of [|, |. Let us analyze M in Eq. fl2j) with S in Eq. Q. 
(a) {i k = 4 = Jk = j'k = 0}. 

Let us take a diagram for W = S*S. Let V\ and v 2 be the vertices inside S, which 
are connected by the propagator 

i r dpo i (3 1} 



V J-oo 2vr P 2 + i0+ 
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(b) {i k = 4 = 1, j k = j{. = 0} and {i_ k = z'_ k = 1, j- k = j\ = 0}. 

We first deal with the case {i k = i' k = 1, j k = j k = 0}. We take out the diagram 
for W = S*S, which is obtained from W above as follows. Remove the propagator 
( p.lD , connect n= i ;k , Eq. (|2.9|) , to the vertex V\ in S, and connect n '=i;k to v 2 - 
Here 0„-i ;k [</> n '=i;k] designates that, in Eq. (|2.9|) , iK k n=1 [iK kn , =1 ] operates on 0„ =1 
[0 n / =1 ]. We pick out from Eq. (|2"Ij|), 

JV™' = N«» = n . (3.2) 

Here and below, we suppress the suffix "k", whenever no confusion arises. In S*, 
Njfl' = N™ = 1. Inserting N$ N*?' = n into Eq. (gj with Eq. (£§), we obtain 

(«) = ^3 1 = ™b(p) • ( 3 - 3 ) 

Here n B {p) = l/(e /3p — 1) is the Bose distribution function and the angular brackets 
denotes the statistical average, 

V°° e>-P n PO 

/O \ = ^n=0 c 

\"n/ - voo « • 

2^n=0 

Then, in A/" in Eq. (|2.4|) , the portion corresponding to Eq. (^) turns out to 



2pV nfl(p) = ^1^^0)2^^)^^), (3.4) 

where l/(2pV) has come from zii^ n '=i^k,n=i in Eq. (|2.9|) with Eq. ( |2.11| ). It is to 
be noted that in K's, Eq. ( |2.11| ), may be dealt with just as in vacuum theory, 

— 1/2 

so that we ignore throughout this paper. 

0-k = i'_ k = !, J-k = jL v = 0}. 

The relative diagram to the above diagram for W = S*S, same as above W except 
that n =i;-k (0n'=i;-k) is connected to the vertex t> 2 (t>i), yields, in place of Eq. (|3.4|), 

i r ^0(-po)27r<y(P 2 )n fl (p). (3.5) 
Adding Eqs. (Q, (gl|), and flPD , we extract 

+ 2nn B ( Pll )6(Pl) = iD n (P k ) 

= ^(PO+zMT^Pk). (3.6) 



P k 2 + i0- 
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Here iDjj and iD^p stand, respectively, for the T-independent part (the first term 
on the left-hand side (LHS)) and the T-dependent part (the second term) of iDn. 

(c) {i k = j k = 0, 4 = j k = 1} and {i_ k = j_ k = 1, i'_ k = f_ k = 0}. 

In order to extract the contribution of {i k = j k = 0, i' k = j k = 1}, we take a 
diagram for W = S*S in Af, Eq. ( |2.4|) , where n '=i;k in S is connected to the vertex 
v\ in 5* and (fi n '=i;ii in S* is connected to the vertex i> 2 in 5*. 

We pick out from Eq. Q2.9| ) and from the form of S 1 *, 

at™/ jv?/ = iv n r +1 ^r +1 = n + 1 . 

Inserting into Eq. (|2.4|) yields 

l + n^l+n B (p). (3.7) 
Then, in M in Eq. ( |2.4j ), the portion under consideration takes the form 

7 /l %o) 27r{1 + 5(p2) • (3 - 8) 

{i_ k = j_ k = 1, z'_ k = = 0}. 

We consider the relative diagram for W = S*S, which is the same as above except 
that n =i;-k in S is connected to the vertex v\ and n =i ; -k in S* is connected to the 
vertex t> 2 . In a similar manner as above, we have 

i r f^(-p ) 27rn B (p) 5(P 2 ) ■ (3.9) 
Adding Eqs. fl3.8|) and ( [Of) , we extract 



2tt [0(p kO ) + n s (p k )] 5(P k 2 ) = iD 21 (P k ) . (3.10) 

(d) {i k = i' k = j k = j' k = 0}, {i k = i k = 0, j k = j k = 1}, and {i_ k = z'_ k = 
0, j- k = j'_ k = l}. 

In a similar manner as in (a) and (b) above, we extract 

+ 2nn B ( Pk )5(Pl) = iD 22 {P k ) (3.11) 



P k 2 - «'0 H 

= {iD xl {P k )y 
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(e) {ik = jk = 1, 4 = 3k = 0} and O-k = J-k = 0, z'_ k = jl k = 1}. 
In a similar manner as in (c) above, we extract 

2n[e(-p ]l0 )+n B (p ]l )] 5(P k 2 ) = iD 12 (P k ) 

= iD 21 {-P k ). (3.12) 

The forms of Dij(P) (i,j = 1,2) defined above are nothing but the thermal 
propagators in the Keldish variant of RTF, which is defined on the time path C, 
— oo — ► +oo — > — oo — > — oo — i/3, in a complex time plane. The above derivation 
shows that the suffix "1" of Dij stands for the vertex in S and the suffix "2" stands 
for the vertex in S*. On the other hand, in RTF, the suffix "1" stands for physical 
or type-1 field and "2" stands for thermal-ghost or type-2 field. 

Let us turn to identify the vertex factors. We take the interaction Lagrangian 
density, 

C int = g^/£\ + \<j>*/e\ . (3.13) 

Then, a (<//) vertex in S receives the factor ig (iX), and then a (<//) vertex 
in S* receives the factor —ig (— iX). This again is in accord with RTF, where ig (—ig) 
and iA (— iX) are the factors which are associated with, in respective order, and 
(//'-vertices of type-1 (type-2) fields. 

Repeating the above procedure, we finally obtain 

(3.14) 

Here A represents the thermal amplitude in the Keldish variant of RTF for the forward 
process, 

E ^( p i) + E ^(Q 3 ) - E <W) + E , 

i=l j=l i=l j=l 

where $i ($2) is a type-1 (type-2) field. The thermal amplitude A is diagrammed in 
Fig. 2. As we have assumed that W = S*S represents the connected diagram (cf. 
above after Eq. ( |2.7|) ), the diagram for A is connected. 
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Each loop momentum P in A accompanies 



In the large V limit the LHS of Eq. ( |3.14| ) becomes 
and Eq. ( |3.15|) becomes 

So far, D in Eq. fl2.5| ) does not participate; T> = 1 (cf. Eq. ( |2.7| )). The role of T> will 
be discussed below. 



3.2 Analysis of mode-overlapping diagrams, ^k + ^k+A+jk > 4 

Above derivation of the thermal-reaction-rate formula is not complete in that we have 
only considered the cases where + z k + jk + j k < 2. When generalized self-energy 
parts are involved in = 5*5", ik + i k + j k + ji > 4. [We call the diagram with 
+ ^k + A + Jk — 4 ^e mode-overlapping diagram.] As mentioned in Sec. I, a 
complete analysis of the classes of diagrams that leads to RTF diagrams including 
thermal propagators with n (> 2) thermal self-energy insertions is still lacking. In this 
subsection, dealing with mode-overlapping diagrams, we shall complete the derivation 
of the thermal-reaction-rate formula. We shall show at the same time that there is 
no finite-volume correction to the formula. 

For illustration of the procedure, we start with analyzing the diagram (for W = 
S*S) with {i k = jk = i'k = Jk = !}• Let us focus our attention on with mode k. 
Both in S and in S*, there are one "absorber vertex" (v[ and t>2 in Fig. 3 below) and 
one "emitter vertex" (^i and v' 2 in Fig. 3). 

From S*S, pick out the factor, 

N nn' N nn> = n 2 ^ 

where and below, the suffix "k" is dropped whenever no confusion arises. In J\f in 
Eq. ( |2.4|) , we have, in place of Eq. ( |3.3|) , 

in 2 ) = 2n 2 B + n B 
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= n B (l + n B ) + n 2 B , (3.17) 

where n B = n B (p). 

The first term on the right-hand side (RHS) of Eq. ( ft.lTQ goes to 

{27r6(p )n B (p) 5(P 2 )}{2n6(p )[l + n B (p)] 5(P 2 )} 

= iD ( £\p)iD£ ) (P), (3.18) 

where D^j 21 (P) = 6(±p ) Di 2 /2i{P)- The (part of) thermal propagator iD[ 2 \p) 
[iDr£\p)] is diagrammed in the double-cut diagram for W = S*S, Fig. 3 (a), as the 
line that connects the emitter vertex v' 2 [i>i] with the absorber vertex v[ \v?\. 
The second term of Eq. ( P-17| ) goes to 

2n6(p )n B (E)5(P 2 )} 2 

= iD^(+\p) lD g ){+ \P). (3.19) 

iD u (P) [iD 22 \P)\ is diagrammed in the double-cut diagram, Fig. 3 (b), as the 
line that connects the emitter vertex v\ [v' 2 ] with the absorber vertex v[ \u?\. Thus, 
with obvious notation, n B {l + n B ) part in Eq. ( |3.17| ) "supplies" (++) portion of 
iD 12 (P)iD 21 (P), Fig. 3 (a), and n| part "supplies" (++) part of iD$P (P) %D$ (P) 
in Fig. 3 (b). 

The ( ) portion of iD\ 2 (P) iD 2 \(P) emerges from W = S*S, which is the same 

as Fig. 3 except that {i_ k = j_k = i'^ = jL^ = 1}- Now, v\ and v 2 [v[ and v 2 ] are 
absorber [emitter] vertices. The (+, — ) portion comes from W = S*S with {i^ = 
i-k = jk = j_ k = 1}. This time, V\ and v[ [v 2 and v 2 ] are absorber [emitter] vertices. 
The ( — h) portion comes from W = S*S with {z' k = i'_ k = j k = ji k = 1}, where the 
absorber [emitter] vertices are v 2 and v 2 \v\ and v[]. Adding all these contributions 
to the contribution ( p.!8| ), we obtain Eq. (|3.18| ) with complete iD 12 (P) iD 21 (P). In a 



similar manner, we can find a set of relative diagrams, which, together with Eq. Q3.19Q , 
yield the complete iDu(P) iD 22 (P). 

All the vertices "t>i", u v[" , u v 2 " , and u v 2 ' (cf. Fig. 3) are not necessarily within 
one connected diagram. There is a diagram as depicted, e.g., in Fig. 4. Figure 4 (a) 
[(b)] contains the factor iD$ iD [ 2 f [iD^ ){+) iD%? W ] in Eq. ( gig ) [Eq. ( gl9|) ]. Let 



us inspect Fig. 4 (a). As stated above after Eq. ( |2.7| ), we are considering the case 
where {A,B} S and {A,B} S * [{A} = {$(P)} and {B} = {$(<?,)}] are involved in 
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one connected part of W — S*S. Then, all $'s are in, e.g., the bottom subdia- 
gram in Fig. 4 (a) (and then also in Fig. 4 (b)) and, in the middle subdiagram, only 
constituent particles 0's of the heat bath participate. iD%y(P) is involved in the 
middle subdiagram, which goes to T>, while iD[ 2 \p) is involved in the bottom sub- 
diagram, which goes to Af con - Thus, Fig. 4 (a) is in Af con D with V ^ 1 in Eq. ( |2.8| ) 
with Eq. ( |2.7D . As a matter of fact, Af con here is obtained from W = S*S with 
{^k = jk = 1, i'k = ik = 0} an d ^ i s obtained from W = SqS (cf. Eq. ( |2.5|) ) with 
{^k = jk = 0, z k = j k = 1}. Thus, Fig. 4 (a) does contribute to 1Z in Eq. ([T^) as 
1Z = Neon, which already appears at lower order of perturbation series. As above, it 
is straightforward to find a set of relative diagrams, which, together with Fig. 4 (a), 
yields the complete iDi 2 (P) iD 2 i(P). Similarly one can find a set of relative diagrams, 
which, together with Fig. 4 (b), yields the complete iDn(P) iD2 2 (P)- 

The relevant part of Fig. 4 (b) and its "relatives" sits in A, Eq. (|3.14 ), as a 



(1, 2) component of a thermal self-energy-inserted propagator. Thus, W = S*S with 
{^k = jk = i'k = jk = 1} together with its "relatives" has turned out to take the 
proper seat in A in Eq. ( |3.14| ). 

It is straightforward to generalize the above argument to a generic diagram for 
W = S*S. Let us focus our attention on a mode k. We analyze M in Eq. ( |2.4| ). Let 
0k be (f) in the mode k. In S in Eq. ( |2.9| ), 0k's in the initial state and i k 0k's in 
the final state participate directly in the reaction. In S*, jk (jk) 0k's in the initial 
(final) state participate directly; ik — i' k = jk — jk = n k — n' k . In S, there are 
(z k ) "absorber vertices" ("emitter vertices") and, in S*, there are jk (j k ) "emitter 
vertices" ("absorber vertices"). [Recall that, in the case of Figs. 3 and 4, v[ and v 2 
are absorber vertices and v% and v' 2 are emitter vertices.] 

We pick out from W = S*S, 

ivrnri ivrnri ^ n \ \ 

jj ' ~ (n-i)\(n-j)\i\i'\j\j>\ 

1 7;U(n + i'-i-k) 3 l[(n-k) , (3.20) 



'• ' ' fc=0 fc=0 



where and below the suffix "k" has been dropped. From the form for S, Eq. (2.S), we 
see that the permutation of <p n i (n' = 1, i k ) and the permutation of <p n (n = 1, i k ) 
give the same diagram, and then i^li^l same diagrams emerge. Then ^k'^k'jk'jk' 
same diagrams emerge for W = S*S, which eliminates the first factor on the RHS of 
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Eq. ( ggUp . In M in Eq. (p^ ), we have, in place of Eq. 



fc=0 fc=0 

Here it is convenient to introduce a generating function of H^j,, 

oo 

f(y, *) = £ V n+l '~ l z n e~ m (x = P P = 

n=0 

In fact, from Eq. Q3.21]) , we obtain 



(3.21) 



TTl, I 



1 d 2 f 



y = z = 1 



/ dy v dxi 

From Eq. (|3~22|) with Eq. ( fOTD , it can be shown that 

i' \ /! (j + i'-A;)! 



(3.22) 



TTl, l 



min(j', j') 

E 

fc=0 



(j>-k)\ 



{n B (x)} 



j+i'-k 



(3.23) 



Since i — i' — j — j', we can readily see that Eq. fl3.23p , is symmetric under 

(i, i') <-» (j, j'). Then, without loss of generality, we assume i> j. 
In Appendix A, we show that 



TTl, l 



J 



min(j,/) 



i'! /! 



(n B y+ k (i + n B y'- k 



(3.24) 



min(j, j') 
k=0 



*{c" 3 ,, k {n B ) k }{cy Kj ,_ k {i + n B y- k } 



(3.25) 



Here n B = n B {p) and 



{i-j+k)\ \ k 

In Eq. ( |3.25| ), the factor C k j may be identified to the number of ways of connecting 
j — k (out of j) emitter vertices in S* to i absorber vertices in S, the factor C° i 
to the number of ways of connecting i — j + k absorber vertices in S to i' emitter 
vertices in S, the factor C® k to the number of ways of connecting k emitter vertices 
in S* to j' absorber vertices in S*, and the factor Cy_ k j,_ k to the number of ways of 
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connecting f — k absorber vertices in S* to %' — {% — j + k) = j' — k emitter vertices 
in 5*. Then, in 7Z in Eqs. (^.3[), we have, in place of Eqs. ( |3.18| ) and fl3.19| ), 



min(i,j') 



j-k 



k=0 



X 



^0 

j+k 



j'-k,j'-k 



{iD£>(p)Y'- h 



(3.26) 



This is just a portion of "right" thermal amplitude in RTF. Just as in the simple 
case, {ik = — i' k — j k = 1}, analyzed above, we can find a set of relative diagrams 
for W = S*S, which, together with Eq. ( |3.26|) , leads to Eq. (|3.26| ) with complete D's. 
Among the diagrams that accompany Eq. ( |3.26j ) with complete .D's, are disconnected 
ones like Fig. 4 (a). Such diagrams belong to Af = N CO nT^ with V ^ 1 (cf. Eq. (|2.8|)), 
and then do contribute to 71 in Eq. (|2~3| ) as 7Z = Af con . Connected diagrams that 
accompany Eq. (|3.26|) with complete D's take the proper seat in A in Eq. (|3.14j) . 

Conversely, for any diagram for A in Eq. ( p. 14 ), through the analysis running 
in the opposite direction, one can identify a set of diagrams for W = S*S. The 
analysis made above is so general that no additional comment is necessary on the 
diagrams that leads to A, Eq. ( |3.14j ), which includes thermal propagator (s) with n 
(> 2) thermal self-energy insertion. 

This completes the derivation of the formula ( |3.14j) for the rate of a generic thermal 
reaction taking place in a heat bath of finite volume. Keeping in mind a suitable 
normalization for incident fluxes of $'s, the formula ( ^.14| ) "smoothly" goes to the 
formula for the infinite-volume (V = oo) system (cf. Eq. ( |3.16 )) in the sense that there 
do not exist extra contributions in Eq. ( |3.14| ) with V < oo, which disappear in the 
limit V —>■ oo. Thus, there is no finite- volume correction to the thermal reaction-rate 
formula ( |3.14j ). 

Here we make a comment on gauge theories. Choosing a physical gauge like 
Coulomb gauge, the gauge boson may be dealt with in a similar manner to the above 
scalar-field case. When we adopt a covariant gauge, Faddeev-Popov (FP) ghost field 
comes on the stage. The first summations in Eqs. ( |2.4| ) and ( |2.5| ) are carried out over 
the modes of physical degrees of freedom. This can be implemented by inserting the 
projection operator V onto the physical space on the left side of p in Eqs. ( p.4|) and 
( p.5| ) and sum is taken over {^k^} for all, unphysical as well as physical, modes ct's. 
As far as the ensemble average of physical quantities like reaction rate are concerned, 
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all the role of V is to make [19[] the antiperiodic boundary condition for FP-ghost field 
the periodic one, (pFp(t — x ) = 4>Fp{t, x ), so that the bare FP-ghost propagator is 
the same in form to the scalar propagator. Keeping this fact in mind, we can deduce 
Eq. (|3.14|) , where A is evaluated using standard gauge- field and FP-ghost thermal 
propagators in the covariant gauge. 



4 The Dirac fermion 



We study the case of Dirac fermion. The expression for 5* in Eq. ( p.9|) with Eqs. ( p,10| ) 
and (|2.11| ) is changed accordingly. Let 



(a = ±) be the number of mode-k 



fermion [anti fermion] with helicity a. The combinatorial factor AT 
changed to 



in Eq. (O) is 



Nf 




(a) (<r)' _(cr)_(cr)' N 
AT n k n k AT n k ™k 
JV ,(a),(a)' JV 7 M 7 M' 



n 

a=± 



where — i^' 



n 




(4.1) 



i^P and — i£ 



-t Y . In Eqs. 



and (12 



the summations on n 



k > 



n. 



7&jj. , and are taken over and 1. We assume 
that the interaction Lagrangian is bilinear in fermion fields, which include fermion 
fields constituting the heat bath and possibly nonthermalized heavy fermion fields, 
the counterpart of $'s in Eq. 



4.1 Analysis of non mode-overlapping diagrams 



We proceed as in Sec. Ill A using the same notation. 



3k 



.(<t)' 
3k 



(a) {i 

In place of Eq. ( |3.1| ) , we have 



l k — 3k _ 3k 



0} (a = ±). 



I 7 



which comes from the following contraction in S (cf. Eq. ( |2.9|) ), 



T 



ijj(xi) ip(xi) ■ ■ ■ ip(x 2 ) ip{x 2 ) 



|0> 



iS F (xi - x 2 )(0|T • ■ • i/)(xi) ■ ■ ■ ip{x 2 ) 



10) 



(4.2) 



(4.3) 
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Here ifjtp's in Eq. (|4.3|) come from the interaction Lagrangian L 



int- 



(+) 



(b) Fermion mode with {i^ = i^p = I, = = 0} (a = ±) and its relative. 
We consider the positive-helicity (er = +) fermion mode with {z k 
1, = j k + ^ =0}. In place of Eqs. ( |3.2|) and (|3.3|) , we have, in respective order, 

N f = n\ 



and 



where 



<n 2 > 



el* + 1 



n F {p) 



n F {p) = l/ief* + 1) 

the Fermi-distribution function and (O n ) ee E^=o e_/3npfi n/ SLo e_/3np - We note 
that the contribution corresponding to Eq. Q4.3| ) above is (cf. Eq. ([2~9|)) 



is 



(0| T [ • • • 1p n '=l(y) ■ ■ ■ 1p{x 2 ) 1p n =l(z) • •• ] |0) 

I I 

= -iS F (y - x 2 ) iS F (x 1 — z)(0\ T ■ ■ ■ ^(a?i) • • ■ ^(ar 2 ) ■ ■ 



10) 



Then, the LHS of Eq. fl3.4|) is replaced by 



2pV 



^-n F (p)w«(F) m (+) (P). 



Adding the contribution from the negative-helicity fermion mode with {z k 



(-) _ „•(-)' 



— 7 Mp)5>M(p)bW(p) 



Adding further the contribution from the antifermion modes with {i_ k 
1, J-k = J-k =0} (a = ±) to Eqs. 



(4.4) 



P k 2 + 



and ( |4.4j ), we extract 
27rn F (p k )<5(P k 2 ) 



= ^n(Pk) = ^r ( ?(Pk) + ^iP(Pk) • (4.5) 
(c) Fermion mode with {i^ = = 0, i^p = — 1} (a = ±) and its relative. 
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In place of Eq. (|3.7|), we have 

1 - n F (p) . 

Then, Eq. ( |3.8| ) is replaced by 

h r T^0(Po)2n{l-n F (p)}5(P 2 )l/>. 

Adding the contribution from the antifermion mode with {i^ k = = 1, z^f k = 
j-l =0} (er = ±), we extract 

2n[6(p ) - n F (p k )}5(P*) / k = tS 21 (P k ) . 

(d) Interchanging the roles of S and S* in (a) and (b) above, we obtain, in place 
of Eq. flPTD , 



27rn F (p k ) 5(P k 2 ; 



P 2 - i0+ 



= iSn{P±) = iS<g>(P k ) + iSg\P u ) . 

(e) Fermion mode with {z k °^ = = 1, i^p = j k a ' = 0} (a = ±) and its relative. 

The relevant statistical factor is n F (p). Let us show that the part under con- 
sideration turns out to iSuiPu)- In place of po > portion of Eq. ( |3.12| ), we have 
2rm F (p k ) 5(P^)/ k which seems to be the po > portion of iSx%{P'k)- However this 
is not the case. Within the resultant reaction-rate formula, which is an amplitude 
in RTF, the above factor 2im F (p\ s ) 5(p k )f > k necessarily appears in association with 
a thermal fermion loop (see below for detail). The thermal fermion loop carries an 
extra minus sign, so that we have, for the portion under consideration, 

iS[t\P k ) = 2ir[-n F (p k )]5(Pl) / k . 

Adding the contribution from the antifermion mode with {z^f k = j^l = 0, iL k = 
j-l =1} (er = ±), we extract 

2n[6(-p ) - n F (p k )} / k 5(P k 2 ) = iS 12 (P k ) . (4.6) 

In the process of deduction, iSji (j, I = 1,2) appears in succession. At the final 
stage, sets of (W) = (S*S) turn out to be thermal amplitudes A's (cf. Eq ( |3.14|) ), 
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which includes thermal loops of the fermion ijj. Out of A's, we take a "standard" 
A s : Each fermion loop contains at most one iS 12. (Note that the number of 1S21 in a 
fermion loop is equal to the number of iS\2.) From A s , we take two fermion loops Li 
and L 2 and let i5 < 2i(P) G L\ and iSi\{Q) € L 2 . iS2i(P) iS2i{Q) comes, with obvious 
notation, from S*S = S*(p, q, ...) S(p, q, ...) = W s , where S is the S'-matrix element 
obtained using Feynman rules (in vacuum theory). The S'-matrix element which is 
related to S(p, q, ...) through exchange p q is — S(q, p, ...), where S(q, p, ...) is 
obtained using Feynman rules. Then, we have 



which brings in an extra minus sign into the corresponding thermal amplitude A. 
Observe here that, through the above replacement of S, L\ and L 2 in A s turns out to 
be an one thermal fermion loop L in A. A thermal fermion loop carries a minus sign. 
Then L\ and L2 in A s carries + = (— ) 2 while L in A carries — . In reducing (W) to 
A, the extra minus sign in Eq. ( |4.7|) eliminates one — , being present in A s , and is left 
with one — , which is interpreted as the minus sign associated with L in A, What we 
have shown is that A is a "right thermal amplitude." 

Repeating the above procedure for "parent" A s 's and "children" A's, as "con- 
structed" above, we can exhaust all A's that contributes to the reaction-rate formula, 
and see that they are "right" thermal amplitudes. 

4.2 Analysis of mode-overlapping diagrams 

Let us turn to analyze the mode-overlapping diagrams. Noting that etc. and 
then also v£ etc. take two values and 1, we shall exhaust all the mode-overlapping 
configurations. 



(a) {i^ = i^p = = =1} (er = ±) and its relatives. 

From Eq. (PD, N f = (n M ) 4 (a = ±), which leads to ((n M ) 4 ) = n F . Through by 



W. 



s 



w 



S*(p, q, ...)S(q, p, ...), 



(4.7) 



now familiar manner, we extract 



n F Y: [{27r6(p )uf\P)u ( f\P)} 




(4.8) 



uf and u\7 \uf and up] in Eq. (|4.8|) are attached to the vertices in S [S*]. 
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The relatives, to be analyzed, of the above configuration are {i 



■ (-a) 
Ji 



■ (-a)' 



1} and {i 



W _ „•(-»)' 



■(cr) _ -(-a)' 



1} (cr = ±). The former yields 
n 2 F J2[{27re(p )u^\p)u^\P)} 

<T=± 

x{2^(p ) M f ) (P)uf(P)}; , 

and the latter yields 

n F (l-n F ) £ [{2ne(p )u^\p)uP(P)} 

<7=± 

x{2^(p )«S <T) (P)^ ( r <T) ( J p)y • 

Adding Eqs. (O) and (O), we obtain 



(4.9) 



(4.10) 



iS^ + \P)).(iS[V i+ \P)) 



j j 



+ n F (l-n F ) [{2n9{p,)uf\P)u^\P)} 

<7=± 

x{2n9(p )u^\P)u^\P)} . 
Adding Eqs. Q4.10D and (|4.11|) , we have 

04? (+) (Pk)) jy 0^ )(+) (Pk))„ 



(4.11) 



(4.12) 



Recalling the fact that % and %' [j and j'] attach to the vertices in S [S*], we see 
that Eq. (|4.12j ) is just a portion of "right" thermal amplitude in RTF. Adding an 
appropriate sets of relative diagrams, we can extract Eq. ( f4.12| ) with complete S"s. 



(b) {z 



(+) _ „•(-) 



jjj. = ' = 1} and its relatives. 



;(-) 



Taking care of the anticommutativity of fermion fields, we extract 



2*e(p ) {u£\P) v£\P) - v£\P) xfc\P)\ 
2n9(p ) {u^(P)u^(P)-u^(P)u ( k \P)} 



(4.13) 



Here t^'s [u/s] are attached to the vertices in S [S*]. Simple manipulation yields 



Eq. (4.13) 



tS[P(P k )). . (iS\?(P k )). 



?(+)/ 



-(iSit\P k )). . (^+ } (P k ) 
19 



12J2 



(4.14) 



Adding appropriate relative diagrams, we can extract Eq. ( |4.14| ) with complete S"s, 
which sits on the "right seat" in thermal amplitude in RTF (cf. Eq. ( ]3.14j )). 



(c) 



;(-) 



,•(+) _ ,•(-) 

Jk — Jk 



1 



■ (a)' 



1} (a = ±) and its relatives. 



We extract 



77.7 



2^(po) {u£\P) u^(P) - v£\p) v£\P)} 
E [2^0)^)7^) 



cr=± 



(4.15) 



where the spinors with suffices and 23 [ji, j'2, and j'3] are attached to the vertices 

in S [S*]. 

We shall show that 



Eq. (4.15)=^(P)-Sj3&(P) 



(4.16) 



where 



«2j2 



?(+)/ 



J 3 '3 



(T)(+) 
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(P) 



+ (<><+> ( p)). (*s!J'(P)).. (i 

+ K'(P)) j (i5i?' + '(P)) j (iS™ + >(P)V 



(4.17) 



We shall prove this by running in the opposite direction, i.e., starting from Eq. (|4.16|) , 
we derive Eq. ( 4.15|) . The first term on the RHS of Eq. ( 4.17 ) consists of two terms, 
the one is proportional to n 2 F and the one is proportional to np. The second and 
third terms are proportional to np. (iS^) (Pfj . . may be written as (cf. Eq. (P~4|)) 

(iS{?(P)) hji = -2nn F (p) 5(P 2 ) £ u^(P)u^(P) . (4.18) 



Other 5"s in Eq. ( |4.17 ) may be expressed similarly. Straightforward but tedious 
manipulation shows that the u n F part" of S 1 ^^— <Sj^}| vanishes. Then, in Eq. (|4.16|) , 
we are left with Ll np part", which turns out to be Eq. (4.15). 
The same comment as above after eq. ( [4.14 ) applies here. 



(d) 



;(-) 



Jk 



(+) 



Jk 



(-) 



;(+)' _ .•(-)' 



Jk 



(+)' 



j'k =1} an d its relatives. 
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We extract 



2n9( Po ) {v£\p) v!r\P) - u { r\P) v£\p)} 
2n9( Po ) {u£\P)u^(P)-u ( k \P)u£\P)}} 
2n9(p ) (P) uH (P) - <7 } (P) < +) (P) } 
2^(po) {^(P) < } (P) - U £{P) uW(P)j 



(4.19) 



As in the above case (c), through straightforward but tedious calculation, we obtain 



Eq. (4.19) 



(P)-S^(P) 



CU«2«3«4 ( p\ _|_ C«2U«3*4 / p\ 



(4.20) 



where 



hhhu \ ) 



is ( 2 l\p) 



iS[t\P)). . (iS[t\Pj). . (iS { 2 t\P)j.. V „ 2J 



J4'4 



+ 
+ 
+ 
+ 
+ 



^SiV {+ \P)).. (iS&(P))..(i 

/lll 3 V ' 1232 V 

^W(P)). . 
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t(T)(+), 



(^) 



l v 



J3J2 V ' J4«4 

(TK+), 



iS[P(P)). . {iSg\P)). . K)W(P)). . 

^(P)Y . (iS£\P)).. (iSg>™{P)). 

/ lo ll \ I 7323 V / 1. 



] (P)). . \ I 2 

/«1«4 



<T)(+, (P)Y . (iS<?< +) (P) 

/ 2124 \ 



J4.71 



i 5 



ii 



7 -<?cn(+) 

'^22 



11 

(T)(+) 



'2 '3 



J3./2 



The same comment as above after Eq. ( f4.14| ) applies here. 

There remains following two configurations to be analyzed; 



(+)' 
k 



■(+)' 
» _ 



1} and its relatives and (f) 



(+)' 



;(-)' 



■(+)' 
Jk 



) {* 



;(-)' 



1, % 



to 



' = 1} (cr = ±) and its relatives. The case (e) [(f)] may be analyzed in a similar 
manner as (b) [(c)] above and the "right combination" of thermal propagators is 
extracted. 

As in the scalar-field case, Sec. Ill B, there appear disconnected Af's; Af = Af con V 
with T> 7^ 1. Such cases are treated in a same manner as in the scalar- field case. 

This completes the analysis of all mode-overlapping configurations. 

Conversely, we take a diagram for A in the reaction-rate formula (cf. Eq. (|3.14|) ). 
The amplitude A contains "vanishing contributions," which should vanish. By this 
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we mean the contributions coming from the configurations, in which at least one of 
(h) .(h)' Ah) Ah)' -Ah) -{h)' Mh) -(h)' (h 



±) is equal to or greater than 2. Let us 
show that such contributions really vanish. Suppose that A contains 



n (^(% 



fc=i 



(4.21) 



where (k = 1,2,3) is the loop momentum (cf. Eq. ( |3.15| )) and the suffix %kjk 
stands for the (ik,jk) element of iSu in the 4x4 Dirac-matrix space. In the loop- 
momentum space, there are "points," where R± = R2 = R3 = R = (ro,r). Adding 
the contributions from the five relative diagrams, we have, in place of Eq. (14.21 ), 



«a 3 n (^(% 

perm k=l 



(4.22) 



where summation is taken over all permutations of (71.72.73)- = +/ — when 

(hhh) is an even/odd permutation of (J1J2J3), which is a reflection of the anticom- 



mutativity of fermion fields. We take the case tq > 0. The "type-1 side" of Eq. ( 4.22 ) 
comes from + = 3, and then ij^ > 2 or i[ ' > 2. Then the contribution 
under consideration should vanish. In order to see that this is really the case, using 
the expression ( OS ), we further extract from Eq. ( [4.22 ) 



E 

perm 



a 



hhh 



n 

/c=i 



o*=± 



K), 



(4.23) 



Again straightforward but tedious manipulation shows that Eq . ( 4.23Q is in fact 
vanishes. In a similar manner, we can show that Eq. ( |4.22j ) with tq < also vanishes. 
We can also see that (iS^(R)j . . lll=2 (^12 (-R)); ,• an d its relatives add up to 

V / Zxjl kjk 

vanish. When product of n (> 4) iSi 2 (R) and/or iS^P(R) appears in A, pick out 
three of them and apply the above argument to show that the contribution vanishes. 

The above analysis applies to all other "vanishing contributions," which include 
nf=i {iS2i(Rk)) with its relatives etc. This completes the proof of absence of "van- 
ishing contributions." 
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5 The rate of reactions between the constituent 
particles of the heat bath 

In the heat bath composed of scalar fields 0's, taking place is the reaction, 



(pi) + ... + 4>(Pm) + heat bath 

-> 0(qi) + ... + 0(q„) + anything . 



(5.1; 



where 0's are the constituent particle of the heat bath. One can easily show that the 
reaction rate takes the form, 



V 



X 



A( P ( 2 ) p(2) nW nW.pW p(i) o (2) O^l 



where A is the RTF amplitude for the forward process, 

01 (Pi) + ... + 01 (P m ) + MQl) + - + MQn) 

-> 2 (Pi) + ... + 2 (P m ) + MQi) + - + MQr 

It is worth noting that Eq. ( |5.2| ) may be rewritten as 



V 



n 



m i 

n 



Li=l 



dp, 



i0 



VJ 2vr 



\V J 2tt 



.4 



.4 



bubble 



(5.2) 



(5.3) 



(5.4) 



The RHS, Abubbie, is a no-leg thermal amplitude, in which no summation is taken over 
Pi (i = 1, ...,m) and q^ (j = 1, ...,n). 

Generalization of the above result to the theories with gauge bosons and/or 
fermions is straightforward. 
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6 Detailed balance 



In this section, on the basis of the generalized reaction-rate formula, Eq. ( |5.2|) , we 
derive the detailed-balance formula through diagrammatic analysis. 

The purpose of this section is to show that the rate ( |5.2|) for the process ( |5.1| ) is 
equal to the rate for the inverse process to ( |5.1| ). [For the case of theories with gauge 
bosons and/or fermions, the same result is obtained.] This is well known for the cases 
of decay- and production-processes, which corresponds to m = 1, n = and m = 0, 
n = 1, respectively, in Eq. ( |5.2| ). 

Take a diagram for A, Eq. ( |5.2| ), and let Ni and N 2 be the number of iD 2 \S and 
iDi2's, respectively, which is involved in A, 

iVi N 2 

II HhiiRj) II iD 12 {R Nl+k ) . (6.1) 

3=1 k=l 

By cutting all the lines iD^s and iD 2 \\ we divide A into one or several "type- 
1 islands" and one or several "type-2 islands". Here, the type-1 (type-2) island is a 
"maximal" amputated subdiagram of A, which consists of only type-1 (type-2) vertices 
and of the propagators iD^s (iD 22 's) connecting them. Then, a type-1 (type-2) 
island includes no type-2 (type-1) vertex. A type-1 (type-2) island is connected by 
iD 2 iS and/or iD^s to type-2 (type-1) island(s). 

Take a type-1 island and we write its contribution (to A) 

Il(Qj 1 , Qj t n Pin •••) Pit) ■ (6-2) 

Here {Pi k , 1 < k < £} is a subset of {Pi, 1 < i < m} on the LHS of Eq. (|Q| ) and 
{Qj k , 1 < k < £'} is a subset of {Qj, 1 < j < n} on the RHS of Eq. ( [5.3|) , where £, 
£' > 0. This type-1 island is connected by S\(> 0) propagators iD 2 is and s 2 (> 0) 
propagators iDi 2 s to one or several type-2 islands. With the help of the identity, 

D 21 ( J R)=e^ D 12 ( J R), (6.3) 

and the momentum-conservation condition, we obtain, for iD's that are attached to 

ft iDniRj) f[iD 12 (R sl+k ) 

j=i k=i 



exp (3 



k=l k=l 



j=i k=i 
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We now take a type-2 island, whose contribution is written as 

-^2 \P' ii j • • • > P%i ) Qji ; • • • > Qjfi , 



(6.5) 



where {Qj k , 1 < k < £'} is a subset of {Qj, 1 < j < n) on the LHS of Eq. ( f5.3| ) and 
1 < fc < is a subset of {P h 1 < z < m} on the RHS of Eq. Q. £ (f ) here 
is not necessarily equal to £ (£') in Eq. ( |6.2|) . In a similar manner as above, in place 
of Eq. (|6.4j), we have, with obvious notation, 



exp [3 



k=l 



fc=l 



fc=l 



fe=i 



For all the islands, we make the above replacements, i.e., the LHS of Eqs. 
and (El]) are replaced with respective RHS. Through this procedure, each iD 2 ± and 
each iD\2 in Eq. (|6.1| ) is "used" twice. Then we obtain 



Eq. 




J=l 3=1 

xl[iD 12 {R j )l[iD 21 (R Nl+k ). 

j=l k=l 

Now we note that the propagators in Xi's (X 2 's) are iD n 's (iD 2 2 , s), and vertices 
in Xi's (X 2 's) are i\ [cf. above after Eq. (|3.13|) 1. Then, using the relation (|3.11|) , 

[iD n (R)Y = tD 22 (R) , 

and [iX]* = —iX, we easily see that| 

■■■1 Qjiii •••> Pit) 

= I 2 (Qj 1 , Qj e ,; P h , PjJ . (6.7) 
§A comment on QCD (QED) is in order. As to the 4-gluon vertex, when compared to the scalar 
theory, no new feature arises. Let Vj (i = 1, 2) be the factor that is associated to a trigluon vertex in 
a type-i island. V, is real and V2 = — Vi. Then, in place of Eq. (6.7), we have X\* = (~) N l2 with N 
the number of trigluon vertices in X\. Since A in Eq. (5.2) contains even number of trigluon vertices, 
Eq. (6.8) holds unchanged. Let us turn to analyze the quark-gluon vertex. In a standard notation, 



the factor associated to a quark-gluon vertex in a type-1/2 island is ±.ig^^T a . Taking trace, in A in 
Eq. ( |5.2| ), of the products of 7-matrices and of color matrices yield a real function of P's and Q's. 
Then, (ig)* = -ig leads to Eq. (E|). To sum up, Eq. @ holds for QCD (QEC). 
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Here we note that, from the first-principle derivation above, it is obvious that, to 
any order of perturbation series, the amplitude A in Eq. ( |5.2p is real, provided that 
all the contributing diagrams are added. This fact, together with Eq. ( p.7|) , shows 
that 

A ( P m p(2) nW pW p(i) o (2) O^l 

/i^-ri ,...,.r m , , ...,t^ n , r x , ...,.r m , ,...,^ n ) 

= exp /3 



i=i i=i 



X 



/UO {2) o( 2 ) p (1) pW-nW o« p (2) p( 2 )l ffi 



)(2) p(l) 



(1) p(2) 



Using Eq. (|6.3|), we obtain 



e Pn n B {Pi) = l + n B ( Pi ), 
e~^{l + n B (qj)} = n B ( qj ). (6.9) 

Substituting Eq. (|6.8|) into Eq. ( |5.2|) and using Eq. ( |6.9|) , we finally obtain 



(6.10) 



Here, the LHS is the rate of the thermal reaction ( |5. 1| ) while the RHS is the rate of 
its inverse process 

0(qi) + ... + 0(q„) + heat bath 

-> 0(pi) + ... + 0(p m ) + anything . 



Equation (|6.10|) is the desired detailed-balance formula. 



7 T —> limit and Cutkosky rules 



In this section, we show that, in the limit, T — > 0, the reaction-rate formula ( p,14| ) 
reduces to the formula that is obtained using the Cutkosky rules. Then, in the case 
of m = 2 and n = 0, Eq. ( 3.14Q goes to the optical theorem and, for m = 2 and n = 1, 
Eq. ( |3.14|) goes to the Mueller formula |20| for inclusive reactions. 

In the previous section, for a given diagram for A in Eq. (|3.14| ), we have defined 
a set of "islands" . The islands in the set may be classified in two groups. The first 
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group consists of the islands, which contains at least one external vertex. Here the 
external vertex is the vertex, in which or from which the external momentum flows. 
The second group consists of the isolated islands, which have no external vertex. 

Let us take the scalar field theory and investigate zero-temperature limit (T — > 
0) of the reaction-rate formula, Eq. ( 3.14] ). [Again, generalization to other theo- 
ries is straightforward.] In this limit, iD 2 \(P) — > 2h6{j)q) S(P 2 ) and iD 12 (P) — > 
2n6(— p ) S(P 2 ). It can readily be seen that, due to momentum conservation, X\ 
and X2, Eqs. (|6.2j ) and (|6.5|) , corresponding to the isolated islands vanish. Then, 
the nonvanishing amplitude A contains only the islands belonging to the first group. 
Thus, we obtain 

s ATi 

a = n[ 2 ^o)^)]n^({^;{Qw 

j=i »=i 

N 2 

xUMiQh; {Ph), (7.i) 

where {P}i etc. denotes the subset of Pi, P m , which flow in the ith "type-1 island" 
etc. {P}i U {Q}i and {Q}j U {P}j are not empty. In Eq. (]7.1|), the direction of all 
the s momenta, i?'s, each of which connects a "type-1 island" and a "type-2 island," 
is taken to flow from the "type-1 island" to the "type-2 island" . As noted before, the 



diagram representing A in Eq. ( 71 ) is connected. 

The RHS of Eq. ( [7.1] ) is just the quantity, which is obtained by applying the 
Cutkosky rules [IT] (in vacuum theory) to the present case. As a special case, consider 
Eq. (|7.1|) with m = 2 and n — 0. Since the particle represented by <fi is stable at 
T = 0, in Eq. Q, N 1 = N 2 = l and {P} t=l = {P} J=1 = {P h P 2 }. Thus Eq. (]71 



is the optical theorem in vacuum theory. Similarly, for m = 2 and n > 1, Eq. ( |7.1| ) is 
just the (generalized) Mueller formula for the inclusive process, 

$(pi) + $(p 2 ) -> $(qi) + ... + $(q n ) + anything. 



8 Thermal cutting rules 

In view of controversy mentioned in Sec. I, we survey in this section the discus- 
sions made in the past for the thermal Cutkosky formula and thermal cutting rules. 
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Although no new result is involved here, it is worth pigeonholing the issue. The 
Cutkosky formula [IT] in vacuum theory is the formula that relates the imaginary 
or absorptive part of an amplitude A to the sum of cut amplitudes Scuts B>( cnt \ For 
simplicity, in this section, we take a self-interacting complex scalar field theory. Gen- 
eralization to other theories are straightforward. i?( cut )'s are constructed from A by 
so cutting the propagators iD's in A that A is divided into As and As*, which are 
amputated. Here As is a part(s) of A and As* is the complex conjugate of the am- 
plitude that is obtained from A by removing As and ID's. Cutting the propagator 
iD{P) makes iD{P) 

2tt6{± Po ) 5{P 2 - m 2 ) , (8.1) 

where the upper (lower) sign is taken when P flows from a vertex in As {As*) to a 
vertex in As* {As). When the Cutkosky formula is applied to a forward amplitude 
A, we see that 1mA is proportional to the corresponding reaction rate, where cutted 
propagators represent the (on-shell) particles in the final state. 

Kobes and Semenoff (KS) |TJ|] were the first who generalized the Cutkosky formula 
to the case of RTF. Namely they obtained the formula that relates the imaginary part 
of a thermal amplitude to the sum of "circled amplitudes," each of which corresponds 
to the "circled" diagram that includes the so-called circled and uncircled vertices. The 
first paper of [12| discusses general thermal amplitudes and the second one discusses 
physical amplitudes, i.e., amplitudes with all external vertices being of type 1. In 
the sequel, unless otherwise stated, we shall restrict our concern to the physical 
amplitudes. The thermal Cutkosky formula deduced in [I2| may be written in terms 
of thermal amplitudes in RTF: 



Im 



iG(P«,..-,pW 
1 



' - 1 n 

~ J2G(P^,...,P^). (8.2) 

ill •••)*n=l 



Here G{P{ , ■ ■ -,Pj^) stands for the (amputated) thermal amplitude with type-i,- 
{j = 1, ...,n) external vertex in which or from which Pj flows. In Eq. (^2), the sum 
J2' stands for taking summation excluding ii = ... = i n = 1 and %\ = ... = i n — 2. 
Note that, as a matter of course, in G, sum is taken over the types (1 and 2) for all 
internal vertices. 
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KS then generalized the notion of cuttings. Comparison of iD 2 i(P), Eq. ( |3.10| ), 
and iD 12 (P) = iD 21 (-P), Eq. ( gTg) , with Eq. (|0|) leads them to regard iD 12 and 
iD 2 i in G"s on the RHS of Eq. (|8.2j ) as cutted propagators. Through cuttings, each G 
is divided into several pieces. KS then introduced a notion of cuttabe and uncuttable 
diagrams. The former diagram is the diagram that does not include isolated island(s) 
(cf. Sec. VII) while the latter diagram includes at least one isolated island. Note that, 
in the case of vacuum theory, all the diagrams are cuttable ones, which motivates KS 
to introduce the above definition. Thus, the terminology "uncuttable" sounded quite 
natural at the time of its introduction. In spite of the fact that this is a matter of 
definition, existence of uncuttable diagrams has aroused controversy. 



Kobes analyzed |TB[] retarded Green functions in terms of circled diagrams. As to 



the usage of "cuttings" , "cuttable" , and "uncuttable," he followed |12 



Jeon analyzed |14| two-point functions in imaginary-time formalism. Continuing 
to the real energies, he discussed thermal cutting rules. His definition of cutting 
is the same as in JT2I], i.e., the propagators iD 12 and iD 21 are regarded as cutted 



propagators. No mention was made on the cuttable and uncuttable diagrams, but no 
doubt that he supposed all diagrams to be cuttable. 

Bedeque, Das, and Naik analyzed |15[ the imaginary part of thermal amplitudes 
(physical and "unphysical) from the same starting formula as in ]I2] , but with different 
route. Recall that the propagator iDjk (j, k = 1,2) connects a type-j vertex with 
a type-k vertex. iDjk is defined to be a cutted propagator if and only if one of the 
type-j and type-A; vertices is of circled and another is of uncircled (cf. the first paper 



of [|T2[). They then showed that the imaginary part of a thermal amplitude is written 
as the sum of cuttable diagrams, in the sense of KS stated above. In each cuttable 
diagram, connected subdiagram(s) at one side of the cut line contains only uncircled 
vertices (external and internal) while connected subdiagram(s) at the other side of 
the cut line contains only circled vertices. As was pointed out in |16| , however, each 
connected part contains in general propagators that are proportional to the on-shell 
factor S(P 2 — m 2 ). Of course, in the zero-temperature limit, their formula as well as 
KS's one reduce to the Cutkosky formula. 

Gelis extensively analyzed [|16j thermal cutting rules for various formulations of 
real-time thermal field theory. As to the usage of "cuttings", "cuttable", and "un- 
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cuttable," he followed [ |12]j . 

Cutting rules for thermal reaction-rate formula are discussed in [^], |3|, 01, ^ [7]. 
Note that, as mentioned above, in vacuum theory, the cutted propagator, Eq. (|8.1|) , 
corresponds to the (on-shell) final-state particle. The thermal cutting rules introduced 
in [Q, |3], [|, [^, |6], [7] is a generalization of this fact. As we have seen above, iG \ 2 (which 
collectively denotes iP>\ 2 and iS\ 2 ) [^21] consists of two parts, the one comes from 
the particle [antiparticle] in the initial state and another comes from the antiparticle 
[particle] in the final state. While iG\ x and iG 22 , the T-dependent parts of iGu 
and iG 22 , come from the interplay of initial-state (anti)particle and the final-state 
(anti)particle. We recall that each of the thermal propagators iGu and iG 22 consists 
of two parts, the T = part iG^ and the T-dependent part iG^ T \ Then, A in 
Eq. ( |3.14|) or ( |5.2|) is divided into 2 N contributions, where N is the number of iGu's 

(T) (T) 

and iG 22 s. Above observation leads us to regard iG\ 2 , iG 2 \, iG n , and %G 22 as the 
cutted propagators. 

Through the applications of the above cutting rules, A is divided into several 
subparts. Each subpart contains only type-1 vertices or only type-2 vertices. The 
former (latter) belongs to S (S*) in (S*S). The cuttings work as follows. The line that 
cut iG\ 2 (P) with pq > (po < 0) is the initial-state particle (final-state antiparticle) 
cut line. The line that cut iG 2 \{P) with p Q > (p < 0) is the final-state particle 
(initial-state antiparticle) cut line. The line that cut iG u (P) [iG 22 (P)\ is the initial- 
state cut line and the final-state cut line in S [S*] and, in S* [S], an one extra spectator 
particle with P is. For the line that cut iGjj(P) with p > (p < 0) is the initial- 
state particle (antiparticle) cut line and the final-state particle (antiparticle) cut line. 
For the cut line on iG 22 (P), similar statement holds. 

It is quite obvious that the "cutting rules" introduced above for thermal reaction 
rates may be used for general thermal amplitudes evaluated in the Keldish variant of 
RTF. 

Finally, it is worth mentioning that it can easily be seen from Eqs. ( p. 14 ) and 
( P-2[ ) that the RHS of Eq. ( |8.2j ), which represents the imaginary part of a physical 
amplitude, is a sum of various reaction rates times corresponding kinematical factors. 
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Appendix A Proof of Equation (3.24) 



Here we prove the identity Eq. ( |3.24j ). We expand the RHS of Eq. ( |3.24| ) in powers 
of 7ig(x)(= £) to obtain 

t Q (i-j + k)\(j'-k)\ y k ) 1 +u 

"fro h (i-j+£)\£\(j-ey.k\(j'-k-£f ' { - } 

where i > j. Comparing Eq. ( |A. 1| ) with Eq. ( |3.23|) , we see that it is sufficient to show 
that 



where 



^ = %; 3 , , (a.2) 



min(i,i'-fc) a\ a\ 

" to £\{i-3+£)\{j'-k-£)\{j-£)\' [ - > 
kQi,i' = U + i'-k)\ 

{i> - k)\ (j' - k)\ ' 1 ' J 



Here we define two functions, 



3 



fc=0 



3 



fc=0 

It can easily be shown that F's and G's satisfy the same differential equation, 

i^(x) = F^Zlix) + j F^ix) , (A.7) 

j-G):' r {x) = GT^rX^+jG^^ix). (A.8) 
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From Eqs. flPQ , flA~!p with Eqs. ( |OD and ([0] ), we obtain 



i^/,(0) = GJ» = ^^ (A.9) 

f; : S( x ) = ( ;%(.,■) ^- ]T . (a. io) 



We see from Eq. ( |A.7| ) [Eq. ( |A.8| )] that F^,(x) \G l -^,(x)\ may be obtained from 



Fy (x) [GY Q (x)) in Eq. (Q) with i' < i', j < j, and K;!.,(0) [0^,(0)] in Eq. Q 



with < i', j < j, j' < j' . Since F's and G's subject to the same set of equations 
fl/TTp - (rA~T0|) , we conclude that 



F-'}(x) = Gl),(x), 



which proves Eq. ( |A.2j ). 

Q.E.D. 
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FIG. 1. Two examples of double-cut diagrams for the transition probability W = 
S*S in vacuum theory. Dashed lines are the final-state cut lines while the dotted 
lines are the initial-state cut lines. The left side of the cut lines represents the 
S- matrix element, S, while the right side does S*. The line that is cutted by the 
final-state (initial-state) cut line represents a particle in the final (initial) state 
in S. The lines cutted by the initial-state [final-state] cut line include those 
corresponding to {A} [{B}] in Eq. ( |2.2|) . The group of lines on top of diagrams 
stands for spectator particles, (a) Both S and S* are connected. In addition 
to the spectator particles mentioned above, additional spectator particles are 
in S*. (b) S is connected while S* is disconnected. Note, however, that S*S is 
connected. 

Fig. 2 Diagrammatic representation of the thermal amplitude A in Eq. (|3.14j) . 

Fig. 3 Double-cut diagrams for W = S*S, which yields (a) iD$ (P) iD$ (P) and 
(b) iDf l ){+ \P)iD t S ){+ \P). Here P = (p,p). 

Fig. 4 Double-cut diagrams for W = S*S, which yields (a) iD$(P)iD£\P) and 
(b) tD^ )i+ \P)iDg H+ \P). Here P = (p,p). 
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